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Summary 

In the absence of specific high-affinity agonists and 
antagonists, it has been difficult to define the target 
genes and biological responses attributable to many 
of the orphan nuclear receptors (ONRs). Indeed, It ap- 
pears that many members of this receptor superfamily 
are not regulated by classical small molecules but 
rather their activity is controlled by interacting cofac- 
tors. Motivated by this finding, we have developed an 
approach to genetically Isolate specific receptor- 
cofactor pairs In cells, allowing us to define the biolog- 
ical responses attributable to each complex. This is 
accomplished by using combinatorial peptide phage 
display to engineer the receptor Interacting domain 
of each cofactor such that it interacts selectively with 
one nuclear receptor. In this study, we describe the 
customization of PGC-1a and its use to study the biol- 
ogy of the estrogen-related receptor a (ERRex) in 
cultured liver cells. 



Introduction 

The ONRs are a subset of the larger nuclear receptor 
superfamily of transcription factors, for which no physi- 
ologically relevant llgands have yet been identified. For 
some ONRs, it has been possible to develop useful syn- 
thetic agonists and antagonists enabling the elucidation 
of their function in vivo. For others, it has been extremely 
difficult to develop small molecule regulators, and con- 
sequently, progress in defining their functional roles 
has been impeded. Some of these problems have 
been mitigated with the advent of siRNA technology, 
wherein knockdown of receptor expression can be 
used to substitute for antagonists. However, the consti- 
tutive activity of many ONRs is so low that even quanti- 
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tative knockdown does not provide the dynamic range 
needed to identify target genes. Thus, there is a clear 
unmet need for approaches with which to positively 
and negatively regulate the activity of ONRs in target 
cells. 

We are interested in defining the physiological roles 
for ERRa, an ONR in the NR3B subfamily. Stmcturally, 
the ERRs are most closely related to the estrogen recep- 
tors (ERs), and not surprisingly, significant crosstalk 
between the signaling pathways regulated by these re- 
ceptors has been observed. However, ERRa also plays 
ER-independent roles as a regulator of oxidative phos- 
phorylation, fatty acid oxidation, and lipid handling 
(Huss et al., 2004; Mootha et al., 2004; Sladek et al., 
1 997). Absent a ligand, however, it has been very difficult 
to evaluate the relative physiological and pathological 
importance of the ERRs in estrogen action and in the 
regulation of metabolism. 

Crystallographic analysis has revealed that ERRa is 
capable of adopting a transcriptionally active conforma- 
tion in the absence of any obvious electron density in 
what would be expected to be its ligand binding pocket 
(Kallen et al., 2004). However, its overexpression in cells 
results in only modest activation of transcription. Robust 
activation is observed when ERRa is coexpressed in 
cells with the coactivator PGC-1a (Huss et al., 2002), 
suggesting that in place of a small molecule ligand, the 
transcriptional activity of this and related ONRs may be 
regulated by cofactor availability. Attempts to develop 
small molecule regulators of ERRa have yielded only 
weak antagonists and inverse agonists whose toxicity 
and receptor crossreactivity limit their use (Willy et al., 
2004). 

One approach to selectively increase the transcrip- 
tional activity of a weakly active receptor is to fuse it to 
the strong transcriptional activator VP1 6. However, this 
approach does not recapitulate the physiological activity 
of the native receptor, as negatively regulated genes will 
likely be turned on by this modified protein. Another 
approach is to activate the receptor with one of its known 
cofactors. However, because most of the coactivators 
that interact with ERRa also interact with multiple recep- 
tors and unrelated transcription factors, it has proven 
difficult, even when complemented with siRNA knock- 
down technology, to study the biology of this receptor 
by simply manipulating coactivator levels in cells. Thus, 
we have developed a methodology that has enabled us 
to selectively and effectively regulate ERRa transcrip- 
tional activity. Specifically, we have used combinatorial 
peptide phage display to engineer the receptor interac- 
tion domain of the coactivator PGC-1 a such that it inter- 
acts with and activates ERRa in a highly selective 
manner. We validated the use of this customized coacti- 
vator by using it to identify ERRa-regulated genes in 
HepG2 cells, a study that both confirmed the key role 
of this receptor in oxidative metabolism and uncovered 
additional pathways in which it is engaged. This technol- 
ogy, we believe, will have broad application, allowing for 
the development of customized coactivators for other 
NRs and transcription factors. 
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Figure 1 . Evaluation of the Receptor Selectivity of the ERRa-Selective Peptides Identified by Phage Display 

The interaction of some ERRa-selective peptides with other nuclecir receptors (NRs) was tested in a mammalian-two-hybrid assay by using 
a VP1 6-receptor construct and a peptide-Gal4DBD fusion on a 5xGal4-luciferase reporter. ERR peptide activity is defined as a percent of highest 
ERR activity: <1 0%; +, 10%-25%; ++, 25%-50%; +++, 50%-75%; ++++, 75%-100%; +++++, >100%; and nd, not determined. ERRa-selective 
peptides with minimal interactions with other NRs are indicated by light or dark shading (light, one to two other NRs; dark, no other NR). Abbre- 
viations: ERR, estrogen receptor-related receptor; ER, estrogen receptor; PR, progesterone receptor; GR, glucocorticoid receptor; AR, andro- 
gen receptor, RAR, retinolc acid receptor; ROR, RAR-related orphan receptor ligand binding domain; TR, thyroid receptor, VDR, vitamin D 
receptor; LXR, liver X receptor; FXR, famesoid X receptor, LRH, liver receptor homolog 1; and PPARy, peroxisome proliferator-activated 
receptor y. 



Results 

Generation of an ERR-Selective Coactivator 
Given that ERRa tethers PGC-1 a to target gene pro- 
moters through either of two LXXLL motifs within the 
coactivator's (CoA) NR-interaction domain (Huss et a!., 
2002), we hypothesized that it may be possible to alter 
PGC-1 a's NR specificity by reengineering its NR-interac- 
tion surface such that it Interacts with ERRa at the ex- 
clusion of other receptors. In this manner, we believed 
that it would be possible to functionally isolate this 
receptor-coactivator pair in target cells and define the 
biological consequences attributable to this complex. 

Previously, we demonstrated that It was possible to 
use combinatorial phage display of peptide libraries 
to identify LXXLL-containing peptides that interact in 
a highly specific manner with the coactivator binding 
pocket of individual NRs (Chang et al,, 1999; Hall et al., 
2000). Thus, using purified recombinant ERRa as bait, 
we screened Ml 3 phage display libraries containing 
10® different random or semirandom peptides 12-19 
amino acids in length. The full details of this procedure 
and the characterization of the peptides identified are 
described elsewhere (Chang et al., 1999; Gaillard et al., 
2006). In brief, we first identified LXXLL peptides that in- 
teracted selectively with ERRa but did not interact with 
either ERa or |3 under the same conditions. Peptides 
that showed a preference for ERR were expressed in 
cells as a fusion with the GAL4DBD, and their ability to 
tether different VP1 6-NR fusions was tested in a mamma- 
lian two-hybrid assay. This secondary screen enabled 
the elimination of those peptides that were (1) capable 
of interacting with NRs commonly expressed with 
ERRa in target cells and (2) receptors known to interact 
with PGC-1 a. The results of this analysis, summarized 
in Figure 1 , indicate that some peptides (L3-09, L3-28, 
and L3-80) are highly selective for the ERRs under the 
conditions tested. 

There are three LXXLL motifs within PGC-1 a 
(Figure 2A, L1-L3); however, only L2 and L3 appear to 
interact with NRs. Although L2 is the dominant site of in- 
teraction for most NRs, we and others have shown that 



either L2 or L3 is sufficient to allow full ERRa transcrip- 
tional activity (Gaillard et al., 2006; Huss et al., 2002; 
Schreiber et al., 2003). Therefore, we replaced both of 
the 19 amino acid regions corresponding to each of the 
L2 and L3 motifs within PGC-1 a with sequences corre- 
sponding to the L3-09 peptide (PGC-1 a 2x9, Figure 2A). 
When assayed on a 3xERE-TATA-lucif erase reporter 
(Figure 2B), PGC-1 a 2x9 stimulates ERRa transcriptional 
activity to the same degree as PGC-1 a. Equivalence of 
PGC-1 a and PGC-1 a 2x9 was also demonstrated in 
experiments using ERRP and ERRy (data not shown). 
ERRa transcriptional activity was unaffected by the 
expression of PGC-1 a L2L3M, a mutant in which the leu- 
cines within the L2 and L3 motifs were changed to ala- 
nines. TTius, within the resolution of these reconstituted 
transcription assays, we conclude that PGC-1 a 2x9 is 
indistinguishable from PGC-1 a as an ERRa coactivator. 

We next tested whether PGC-1 a 2x9 could activate 
other NRs, in particular those known to be coexpressed 
with ERRa and/or subject to coactivation by PGC-1 a. 
To this end, we compared the coactivator activity of 
PGC-1 a or the mutants on individual receptors and cog- 
nate reporters In the presence or absence of a ligand as 
required. Of the receptors tested that are known to inter- 
act with PGC-1 a, only HNF4a activity was enhanced to 
any significant degree by PGC-1 a 2x9, resulting from 
crossreactivity of the L9 peptide, whereas the transcrip- 
tional activity of other receptors such as ERa, PPARy, 
GR, LRH-1 , or RXRa was unaffected (Figure 2B). Thus, 
by engineering the LXXLL motifs within PGC-1 a, we 
have been able to develop a highly selective coactivator 
for ERRa. 

Identification of ERRa Target Genes in HepG2 Cells 
by Using an ERR-Selective PGC-1 a 
A role for the ERRa/PGC-1 a complex has been shown in 
the regulation of mitochondrial biogenesis, oxidative 
phosphorylation, and fatty acid oxidation in cardiac 
muscle and bone precursor cells (Huss et al., 2004; 
Schreiber etal., 2004). PGC-1 a, HNF4a, and other recep- 
tors have a well-established role in the regulation of glu- 
coneogenesis in liver cells (Puigserver and Spiegelman, 
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Figure 2. ERR-Setective PGG-1a Preferen- 
tially Coactivates ERR 

(A) Schematic of PGC-1a LXXLL region 
mutants. 

(B) HeLa cells were transfected with the indi- 
cated receptor, cofactor, and luciferase re- 
porter conteuning the indicated response ele- 
ment. One-hundred nanomolar 17p-estradiol 
(E2), 100 nM dexamethasone, and 100 nM 
9-c/s-retinoic acid were used to activate ERa, 
OR, and RXRa, respectively; no hormone indi- 
cates an equivalent volume of ethanol. Re- 
sults are expressed as normalized luciferase 
activity ± standard error of the mean (SEM) 
per triplicate sample of cells. 



2003) , but the role of ERRa in mediating energy metabo- 
lism in these cells has not been defined. As an initial step 
toward addressing this issue, we elected to define the 
target genes expressed in the hepatocellular carcinoma 
cell line (HepG2) that were subject to regulation by the 
ERRo/PGC-l a complex. 

We first determined whether PGC-1 a 2x9 could acti- 
vate physiologic targets when transduced into target 
cells. PGC-la facilitates transcriptional auto upregula- 
tion of the mRNA encoding ERRa (Laganiere et al., 

2004) . Thus, we measured the expression of ERRa in 
HepG2 cells overexpressing either PGC-1 a or the 2x9 
variant. As expected, ERRa mRNA and protein levels 
were strongly induced by both the PGC-1 a and PGC-la 
2x9, but not by either of the control viruses as measured 
by quantitative PCR (qPCR) and immunoblot, respec- 
tively (Figure 3A). We observed no induction of endoge- 
nous PGC-la by either the PGC-la- or PGC-la 2x9- 
expressing adenoviruses (Figure SI in the Supplemental 
Data available with this article online) and proceeded to 
use the targeted PGC-la to define the ERRa transcrip- 
tome in HepG2 cells. 



To identify target genes induced by PGC-1 a through 
ERRa, PGC-1 a, the 2x9 variant, and the controls were 
independently transduced in HepG2 cells and the corre- 
sponding gene expression profiles were analyzed by 
using the GeneChip technology (Affymetrix, Santa Clara, 
CA). The resulting hierarchical clustering diagram (Fig- 
ure S2) illustrates that the vast majority of genes are sim- 
ilarty regulated by PGC-1 a and PGC-la 2x9, differing 
only by degree of induction or repression. Of the genes 
significantly induced, 94% were induced by both PGC- 
1 a and PGC-1 a 2x9, whereas 76% of the downregulated 
genes were repressed by both (Figure SB). This high 
percentage of similarity suggests that PGC-1 a 2x9 reca- 
pitulates the activating capacities of PGC-1 a and that 
a significant number of PGC-1 a target genes in HepG2 
cells are regulated by ERRa and/or HNF4a. However, 
quantitative differences exist between the magnitude 
of induction or repression by PGC-1 a 2x9 compared to 
PGC-1 a. PGC-1 a induces 61 % of the upregulated genes 
more effectively than PGC-1 a 2x9, suggesting that opti- 
mal expression of these genes relies on the additional 
interactions by PGC-1 a, in which PGC-la 2x9 is unable 
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Figure 3. Regulation of Transcription by WT 
PGC-1a and PGC-1a 2x9 in HepG2 Cells 

(A) Induction of ERRa mRNA and protein by 
PGC-la constructs. Enror bars indicate SEM 
of three biological replicates. 

(B) Qualitative summary of genes significantly 
regulated by PGC-la or 2x9 or both. 

(C) Quantitative differences In gene regula- 
tion by PGC-1 a or 2x9 based on hierearchical 
clustering. Expression levels from clusters 
compared between treatment groups using 
unpaired Student's t test (p < 0.002, where 
degree of Induction or repression was differ- 
ent, p > 0.1 96 for similarly regulated clusters). 

(D) Gene ontology (GO) categories that are 
most prominently Influenced by PGC-1 a and 
PGC-la 2x9 related to energy metabolism. 
P values indicating the enrichment of these 
GO terms among the regulated genes were 
conrected for multiple comparisons. 
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to participate (Figure 3C). However, 39% of the upregu- 
lated genes are induced by PGC-1 a 2x9 as well or better 
than PGC-1 a, suggesting that these genes are primarily 
regulated by either ERRa or HNF4a, For 67% of the 
downregulated genes, PGC-la 2x9 acts as a stronger 
corepressor. This may result from the fact that the effec- 
tive pool of PGC-1 a available to ERRa is greater, as other 
transcription factors cannot effectively compete for 
binding to the modified coactivator. 

Analysis of the gene ontology (GO) temis associated 
with the differentially regulated transcripts revealed an 
enrichment of terms related to energy metabolism, In- 
cluding oxidative phosphorylation, tricarboxylic acid 
cycle, fatty acid beta-oxidation, lipid metabolism, lipid 
transport, and electron transport (Figure 3D). Most of 
these pathways have been shown previously to be 
targets of ERRa (Huss et al., 2004; Mootha et al., 2004; 



Sladek et al., 1997). Other categories such as mitochon- 
drial transport, regulation of transcription, and amino 
acid metabolism did not show statistical significance 
in number, but the biological significance of the individ- 
ual genes in each class should not be discounted 
(Figure S3). 

The PGC-1 o/ERRa Complex Is a Key Regulator 
of Rate-Limiting Enzymes Involved 
In Energy Metabolism 

Potential ERRa targets identified in the array described 
above were validated In independent experiments by us- 
ing qPCR after infection with a virus expressing PGC-1 a. 
To differentiate between ERRa- and HNF4a-dependent 
targets, we compared the relative levels of mRNA of 
PGC-1 a-lnduced genes in cells treated with an adenovi- 
rus expressing either a control (scrambled) or ERRa 
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Figure 4. PGC-1 a-Mediated Induction of Key Enzymes of Energy Metabolism Is Dependent on the Expression of ERRa 

(A) Induction of ERRa mRNA and protein by PGC-1 a or PGC-1 a 2x9, and repression by the siRNA to ERRa. 

(B) Genes not dependent on ERRa. 

(C) ERRa-dependent genes. 

Error bars in (A)-(C) represent SEM of three biological replicates. 

(D) Fatty acid oxidation in HepG2 cells, normalized to protein content, with error bars indicating SEM of three independent experiments. 



siRNA. The siRNA effectively inhibits the induction of 
ERRa mRNA and protein by PGC-1 a as compared to 
a scrambled control siRNA (Figure 4A]. Under these 
conditions, it was possible to distinguish known HNF4a- 
activated genes, such as phosphoenolpyruvate car- 
boxykinase 1 (PCK1), cholesterol 7-alpha-hydroxylase 
{CYP7A1), and glucose-6-phosphatase {G6PD) (Fig- 
ure 4B and data not shown) from key genes of energy 
metabolism found to be dependent on ERRa, includ- 
ing isocltrate dehydrogenase 3A {IDH3A) and oxalogluta- 
rate dehydrogenase (OGDH) of the tricarboxylic acid 
cycle, peroxisome proliferator-activated receptor alpha 
(PPARa) and carnitine palmitoyltransferase 1A (CPU A) 
involved in regulating fatty acid oxidation, ATP synthase 
P {fiJPsynp, cytochrome c (Cyt c), and cytochrome c ox- 
idase 4 (COX4) of oxidative phosphorylation and adenine 
nucleotide translocator 1 (ANT1) of mitochondrial trans- 
port. It is important to note that of all the genes we have 
studied thus far to be induced by PGC-1 a 2x9 In HepG2 
cells, PCK7, CYP7A1, and G6PD were the only genes 
shown to be dependent on a factor other than ERRa. 

To validate this approach in other cell types, we tested 
induction of these ERRa target genes in AGS stomach 



cancer cells and U251 glioma cells and found that these 
genes were similarly induced by both PGC-1 a and PGC- 
la 2x9 in an ERRa-dependent manner (Figure S4 and 
data not shown). Additionally, we tested the robustness 
of the customizing process by exchanging a different 
ERRa-selective peptide, L28, for the two L9 peptides 
and observed that the ERRa target genes tested were in- 
duced by the PGC-1 a 2x28 as well (Figure S5). Finally, to 
determine whether the PGC-1 a 2x9 construct is func- 
tionally indistinguishable from PGC-1 a in inducing one 
of the expected pathways, we measured fatty acid oxi- 
dation in HepG2 cells. We found PGC-1 a 2x9 Increases 
oxidation of oleic acid compared to cells expressing 
3-gal or PGC-1 a L2L3M and recapitulates PGC-1 a 
activity through ERRa (Figure 4D). 

Discussion 

It is unlikely that classical small molecule ligands will be 
found for all ONRs. Rather, it appears that coactivator 
availability or posttranslational modifications, such as 
phosphorylation, will emerge as primary mechanisms 
by which the transcriptional activity of some ONRs will 
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be regulated (Hermanson et at., 2002). ERRa appears to 
interact primarily with members of the PGC-1 and/or 
p160 classes of transcription factors (Lu et al., 2001; 
Schreiber et al., 2003). Because these cofactors interact 
with several NRs and with other unrelated transcription 
factors, it has been exceedingly difficult to define the 
biological responses in cells that are attributable to indi- 
vidual NR-CoA complexes. With this problem in mind, 
we developed an approach that has allowed us to genet- 
ically isolate ERRo/PGC-la within cells and identify 
target genes that were regulated by this specific NR- 
CoA complex. This was accomplished by engineering 
the LXXLL NR-interacting motifs within PGC-1 a such 
that they would interact with ERRa at the exclusion of 
other NRs. Although the LXXLL motifs within the known 
coactivators are relatively promiscuous, we were able to 
identify peptides that were highly selective for ERRa by 
screening large peptide libraries in which the amino 
acids flanking the core motif were randomized. We be- 
lieve that this general approach to "customize" coacti- 
vators will be applicable to other NRs, as we have found 
it to be relatively easy to generate LXXLL-based pep- 
tides that display a high degree of receptor selectivity. 
Although it is unlikely that absolute specificity of the 
coactivator for the targeted receptor can be achieved 
with this technology, this problem also exists for natural 
ligands of the NRs. For instance, progesterone is an 
effective ligand not only of the progesterone receptor 
but also of the mineralocortlcoid and glucocorticoid 
receptors. Thus, confirmation of an observed biological 
response to a small molecule "classic" agonist or a cus- 
tomized coactivator will always be required. To this end, 
we have used siRNA in this study to confirm the involve- 
ment of ERRa in the highlighted processes. Through the 
use of these two technologies, comparable to the use of 
classical agonists and antagonists, we believe that it will 
be possible to probe the biology of most any transcrip- 
tion factor-cofactor pair. 



Experimental Procedures 
Plasmids 

Plasmids for pcDNA3-ERRa (Zuercher et al., 2005); VP16-ERa, 
VP16-ERP, VP16-RARa, VP16-RXRa, pRST7-ERa, 3xERE-TATA-lu- 
ciferase, and 5xGal4-luciferase (Chang et al., 1999); LRH-1 (Safi 
et a!., 2005); and VP1 6-LXR and VP1 6-FXR (Hall et al., 2000) were de- 
scribed previously. cDNA encoding ERRa and PPARy were excised 
from pcDNA3-ERRa and SGS-PPARy, respectively, and ligated into 
the VP16 vector (Clontech), The following were gifts: SGS-PPARy 
and SxPPARE-lucif erase (T. Willson), SHP-luc if erase (S. Kliewer), 
DR1 -luciferase and RSV-RXRa (R. Heyman), HIV A1 -4-luclferase 
and pcDNA3-HNF4a (F. Sladek), Flag-PGC-1a (B. Spiegelman), 
VP16-VDR (W. Pike), VP16-TRP (D. Moore), VP16-PRA (D. Wen), 
VP16-GR (J. Miner), VP16-AR (K. Marschke), and VP1 6-RORa-LBD 
(A. Means). 



Generation of the PGC-1a 2x9 and L2L3M Mutants 
PGC-1 a 2x9 was generated by excising the L2 and L3 motifs of PGC- 
1 a and ligating in their place the sequence of the L3-09 peptide 
(PSEDRGELWRLLSVTERQN). This resulted in the replacement of 
the L2 and L3 motifs with an ERR-selective peptide and the insertion 
of a total of seven amino acids conresponding to restriction enzyme 
sites. The PGC-1 a L^L3M mutant, in which the leucines were 
mutated to alanines, was generated by using the QuikChange II 
mutagenesis kit (Stratagene). All mutants were sequenced to ensure 
the fidelity of the resulting constructs. 



Generation of Adenoviruses 

Adenoviruses used to overexpress PGG-1 a, PGC-1 a 2x9, or PGC-1 a 
L2L3M were generated by using the ViraPower Adenoviral Expres- 
sion System (Invitrogen), purified with Adeno-x (BD Biosciences), 
and concentrated with a Centricon YM-50 Unit (Millipore). 

The ERRa and control siRNA adenovimses were generated as 
described above with the following changes: the ERRa siRNA oligo- 
nucleotides (described in Schreiber et al. [2003]) and control siRNA 
oligonucleotides (scrambled sequence) were ligated into pSuper 
(Oligoengine) and the H1-RNA promoter and the target sequence 
were inserted into the adenoviral backbone. 

Mammalian Cell Culture, Transfections, 
and Adenovirus Transduction 

Culture, transfection, and luciferase assays using HeLa (human cer- 
vical carcinoma) and HepG2 (hepatoma) cells were described previ- 
ously (Chang et al., 1 999). For assays Involving homnone receptors, 
cells were treated with vehicle or the following concentrations of its 
respective hormone: 100 nM 17p-estradiol (ERa, ERP), 100 nM 
progesterone (PR- A), 1 00 nM dexamethasone (GR), 1 i^M 5a-dihy- 
drotestosterone (AR), 100 nM 9-c/5-retinoic acid (RARa, RXRa), 
100 nM triiodothyronine (TRP). 100 nM 1 ,25-dihydroxyvitamin 03 
(VDR), 10 ^M 22R-hydroxycholesterol (LXR), and 50 ^M chenodeox- 
ycholic acid (FXR). For transduction of protein using adenoviruses, 
HepG2 cells were infected at a multiplicity of infection (MOI) of 
10-50 for 48-72 hr. In experiments with siRNA and PGC-1 a infec- 
tions, cells were infected with siRNA adenovirus 48 hr prior to 
a second infection with PGC-1 a (or mutants). 

Protein Immunoblotting 

Whole-cell protein extracts were separated on an 8% SDS-PAGE 
and blotted onto PVDF (Millipore). PGC-1 a protein was detected 
with a rabbit polyclonal PGC-1 a antibody (Santa Cruz Biotechnol- 
ogy). ERRa protein was detected by using a monoclonal antibody 
(Gaiiiard et al., 2006). 

Quantitative PCR 

Total RNA was isolated by using the RNeasy kit with RNase-free 
DNase (Qiagen). One microgram of RNA was reverse transcribed 
by using the BioRad iScript cDNA synthesis kit. qPCR reactions 
(BioRad iCycler) were performed with 0.1 ^1 of cDNA, 1 0 ^M specific 
primers, and IQ SYBRGreen supermix (BioRad). The sequences of 
the primers are listed in Table SI . 

Fatty Acid Oxidation 

HepG2 cells were infected at MOt 20 for 72 hr then incubated 3 hr at 
37°C with 1 mM carnitine, 12.5 mM HEPES, 0.5% BSA. 500 ^M so- 
dium oleate, and 1 .0 ^Ci/mi [1 -^'^0] oleic acid (MP Biomedicals). Ox- 
idation end product ^^C-COz released was measured as described 
previously (Koves et al., 2005). 

Microarray Statistical and Gene Ontologic Analyses 
Microarray gene profiling experiments were performed with the 
Human Genome U133 Plus 2.0 Anray chips (Affymetrix). RNA was 
collected 24 hr after Infection. Target preparation, hybridization to 
the Affymetrix HG-U133 plus 2.0 arrays, and scanning were per- 
formed by the Duke Microarray Center. The entire experiment was 
repeated three times. The probe intensities were extracted from 
CEL files by utilizing the Affymetrix Input Engine from the SAS Micro- 
array Solution software (SAS Institute Inc., Gary, NC). After normal- 
ization, the mixed model was applied by running Mixed Model Anal- 
ysis from the SAS Microan-ay Solution (Chu et al., 2002, 2004). 
Significant genes were determined by conducting t tests based on 
the estimated parameters from this model. The specific tests con- 
sisted of differences between p-gal control and three PGC-a variants 
treatment groups (WT, 2x9, and L2L3M). A Bonferroni correction 
was applied across all tests to control the probability of false posi- 
tives to be <0.05. Enrichment in GO terms was detected by using 
the Onto-Express application (Draghici et al., 2003). 

Supplemental Data 

Supplemental Data include five figures and one table and can be 
found with this article online at http://www.molecule.org/cgi/ 
contenVf ull/24/5/797/DCI /. 
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Accession Numbers 

Microarray .eel and .chp files can be found on the NURSA micro- 
Einray database (http://www.nursa.org/datasets.cf m?doi+1 0.621/ 
datasets.02004) and the Gene Expression Omnibus (http://www. 
ncbi.nlm.nih.gov/geo/), accession number GSE5968. 



